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Abstract 


An  analysis  of  a  twisted  turbine  blade  was  performed  to  determine  the  regions  of  maximum 
stress  and  moment  which  occur  on  a  typical  gas  turbine  engine  under  variable  revolutions  per 
second.  An  extension  of  the  Holzer  and  Myklestad  methods  was  developed  to  simulate  shear 
and  moment  effects  occurring  from  the  root  of  the  blade  to  its  outermost  extremity.  Under 
numerical  simulations,  shear  force  and  moments  of  inertia  were  realized  at  changing  revolutions 
per  minute  by  solving  the  Den  Hartog  twisted  blade  equations. 
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1,  Introduction 


The  usual  method  for  determining  the  natural  frequencies  or  critical  speeds  of  shafts  or  beams 
in  bending  is  the  iteration  method  of  Stodola  (Anderson  1967;  Crandal  1956;  Den  Hartog  1956; 
Harris  and  Crede  1961;  Kelly  and  Richman  1971;  Meirovitch  1967;  Timoshenko,  Young,  and 
Weaver  1967;  Tong  1960;  Vernon  1967;  Andronow  and  Chaiken  1949;  Cunningham  1958; 
Minorisky  1962;  Stoker  1950)  either  in  graph  or  numerical  form.  In  addition,  another  extension  was 
introduced  called  the  Holzer  method,  which  has  been  applied  to  flexural  vibration.  In  the  torsional 
problem,  the  angle  and  twisting  moments  are  significant  calculational  quantities;  in  the  flexural 
problem,  there  are  usually  the  deflection,  the  slope,  the  bending  moment,  and  the  shear  force,  which 
are  pertinent  and  relevant  quantities.  It  is  necessary  to  find  the  relations  among  these  quantities  from 
one  section  (station)  to  another  along  a  turbine  blade  of  constant  cross  section  (nontwisted). 

Field  and  point  transfer  matrices  may  also  be  used  to  obtain  relations  that  govern  the  flexural 
motion  and  vibrations  of  turbine  blades.  The  transfer  matrices  are  suitable  for  solution  by  the 
Myklestad  method  (Tse  1978),  which  is  essentially  a  Holzer  procedure  where  the  primary  difference 
lies  in  the  use  of  the  equations  for  bending  and  the  manipulation  of  these  equations  to  solve  readily 
in  terms  of  the  boundary  conditions  for  the  blade  itself;  in  this  particular  case,  a  remainder  slope  is 
used  as  the  test  for  the  selection  of  the  correct  frequency. 

However,  the  Myklestad,  and/or  Holzer  method,  can  be  extended  to  more  complicated  cases,  for 
example,  to  a  twisted  turbine  blade  in  which  the  principal  axes  of  bending  stiffness  turn  through  an 
angle  along  the  blade  length.  Then,  a  vibration  in  one  plane  is  no  longer  possible  (i.e.,  the  motion 
in  the  axial  direction  is  coupled  to  that  in  the  tangential  direction).  Hence,  the  purpose  of  this  report 
is  to  analyze  the  shear  and  moments  of  a  typical  twisted  turbine  blade  as  extracted  from  a  gas  turbine 
engine.  The  information  yielded  in  this  analysis  will  aid  and  assist  the  understanding  of  the 
mechanical  vibrations  anamolies  associated  with  the  Turbine  Engine  Diagnostics  (TED)  project 
(Helfrnan,  Dumer,  and  Hanratty  1995)  to  develop  acoustical  signatures  necessary  for  performance 
assessment  and  prediction. 
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2.  Analysis 


Let  the  axial  direction  be  x,  the  tangential  direction  be  y,  and  the  radial  direction  along  the  blade 
be  z.  First  calculate  the  mass  per  unit  length,  p,,  the  bending  stiffnesses,  and  ly,  in  the  axial  and 
tangential  directions,  and  the  product  of  inertia,  I^^y,  or  stiffness  coupling.  All  of  these  are  variable 
with  the  length  z.  Then  cut  the  blade  into  a  number  of  sections;  take  the  average  value  of  pi,  1^,  ly, 
and  L„  of  each  section  and  consider  these  values  constant  along  each  section.  The  vibration  takes 
place  in  the  x  and  y  coordinates  simultaneously  so  that  there  will  be  eight  equations.  There  is  cross 
coupling  on  account  of  the  product  of  inertia  I,ty  term.  The  eight  equations  are  given  below  without 
derivation  (Den  Hartog  1956). 
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Here, 

S,j,  Sy  represent  the  shear  forces  in  the  x  and  y  directions,  respectively; 

M,t,  My  represent  the  moments  in  the  x  and  y  directions,  respectively; 
m^  is  the  mass  of  nth  section  or  station  of  the  twisted  blade; 

L  is  the  total  length  of  the  twisted  blade; 

I,;,  ly  represent  the  moment  of  inertia  in  the  x  and  y  directions,  respectively; 
and  I,(y  is  the  product  of  inertia. 

The  boundary  conditions  are  as  follows: 

(1)  At  station  zero  (root  of  blade): 

x  =  y  =  x'  =  y'  =  0.  (9) 

The  bending  moments  and  shear  forces  at  the  root  are  unknown.  Take  =  1 ,  My  =  My^,  =  S,,,, 
and  Sy  =  Syo  at  the  root. 

(2)  At  the  free  end,  we  have 
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M,  =  My  =  S,  =  Sy-0. 


(10) 


The  solution  methodology  will  involve  casting  the  set  of  eight  equations  in  matrix  notation  to 
effectuate  a  solution  for  given  values  of  the  frequency,  w.  Tlie  matrix  system  to  be  solved  is  shown 
in  Table  1  in  which  the  coefficients  Cl, ... ,  C16  can  be  defined  as: 


C,  =  -L%/[E(I,Iy-4)] 

(11) 

C2  =  -2CiL/3 

(12) 
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C,3  =  C5/L 

(23) 

Cj4  =  -Cj3L /  2 

(24) 

C,5  =  C3/L 

(25) 

Cig=  -C15L/  2 

(26) 

such  that  the  index,  n,  has  been  omitted  for  convenience. 

The  solution  was  achieved  through  a  matrix  inversion  technique  in  the  IMSL  suite  of  software 
residing  on  the  Cray  C90  in  the  Waterways  Experimental  Station.  Several  values  of  the  angular 
velocity,  w  were  varied  such  as  1,000-9,000  rpm  to  obtain  a  realistic  range  of  operating  conditions 
in  a  gas  turbine  engine. 

3.  Results 

The  distribution  of  the  x-component  of  the  shear  stress  along  various  stations  on  the  blade  itself 
at  numerous  revolutions  per  minute  is  shown  in  Figure  1 .  Here,  the  maximum  shear  takes  place  next 
to  the  root  of  the  blade  and  then  diminishes  toward  the  end  or  tip  of  the  blade  itself,  which  is  very 
similar  to  a  typical  cantelever  beam  under  a  uniformly  distributed  load.  Figure  2  depicts  the 
distribution  of  the  y-component  of  the  shear  stress  at  stations  along  the  blade  for  variable  revolutions 
per  minute  where  again  the  maximum  shear  occurs  immediately  adjacent  to  the  root  of  the  blade  as 
intuitively  expected. 

Figure  3  illustrates  the  axial  distribution  of  the  x-component  moment  of  inertia  for  the  blade  at 
different  stations  and  at  variable  revolutions  per  minute;  clearly,  it  is  seen  that  the  moment  becomes 
a  maximum  at  around  one  third  the  distance  from  the  root  of  the  blade  for  the  largest  revolutions  per 
minute  (or  for  increasing  revolutions  per  minute).  Similarly,  Figure  4  shows  the  maximum  moment 
of  inertia  (y-component)  occurring  near  the  immediate  vicinity  of  the  root  of  the  blade  for  large 
revolutions  per  minute. 
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City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


DEPARTMENT  OF  THE  ARMY 


OFFICIAL  BUSINESS 


I 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  0001,APG,MD 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


DIRECTOR 

US  ARMY  RESEARCH  LABORATORY 
ATTN  AMSRL  IS  Cl 

ABERDEEN  PROVING  GROUND  MD  21005-5067 


